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ABSTRACT 



Much analysis has been done to date on the deformation of helical springs 
due to normal loading. The aim of this study is to design a helical spring that will 
deform under eccentric loading a desired amount due to a given force. Under the 
assumptions of linear stress strain relationships, the spring will be designed in 
terms of its material properties and its geometry. The deformation of the spring 
will be made possible utilizing a Shape Memory Alloy (SMA) active element that 
undergoes phase transformation upon heating above a certain temperature. Two 
models for spring deformation have been considered. In the first model we study 
the differential compression of a spring using SMA wire actuators, and in the 
second model we investigate the bending of an SMA rod placed inside the spring. 
Our efforts were a first step towards the development of a structural skeleton for 
a minimally invasive surgical manipulator. 
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I . INTRODUCTION 

The study of the deformation of helical springs has most commonly been 
limited to those cases due to normal loading [Ref. 1 , 4 , 5 , 6 , 7 , 8 ] . In 
normal loading, a force is applied through the center of a spring, and the 
displacement of the spring is expressed as a function of the load and the 
parameters of the spring involving material properties and geometry. This 
thesis investigates the deformation of helical springs due to eccentric 
loading, with the aim of designing a helical spring that will deform to a 
given shape through the application of known forces. The spring design 
will include the material selection and the selection of helical spring 
geometry in terms of its length, the diameters of the spring and the 
spring coil, and the number of coils. 

Chapter II provides the background needed to study the deformation of 
the two helical spring models, such as the development of Cast igl iano ' s 
theorem, and a brief introduction of Shape Memory Alloy's (SMA's). The 
descriptions of the two kinematic models that have been chosen for this 
study are also included in this chapter. Chapter III provides the 
analysis for spring deformation based on the first model, and Chapter IV 
provides the same for the second model . Chapter V contains computer 
programs that will be used for simulation based design of the proper 
helical spring. 

The motivation of this study is to design a helical spring that will 
be used as the primary structural element of a robotic manipulator for 
minimally invasive surgical applications. In viewing the spring as a 
skeletal element, it is necessary to analyze the deformation of the 
skeleton under the action of external forces, so as to control its 
deformation or motion. The spring deformation will be produced through 
the use of Shape Memory Alloy (SMA) active elements which undergo a phase 
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transformation upon heating above a certain temperature. When these SMA's 
undergo their phase transformation, they change their shape to a 
predetermined form. During this process they will deform the spring they 
are acting on. By properly selecting the arrangement of the SMA active 
elements about the spring, and controlling their phase transformation, a 
properly designed spring can be made to deform to a desired shape. 
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I I . PRELIMINARIES 



Certain preliminary topics need to be reviewed before an analysis of 
the deformation of helical springs due to eccentric loading can be carried 
out. This chapter includes the development of Cast igl iano ' s theorem, an 
introductory examination of the properties of Shape Memory Alloy's, and 
descriptions of the spring models we have chosen for analysis. 

A . CASTIGLIANO ' S THEOREM 

Castigliano' s theorem is a very useful mathematical tool that can be 
used to study the deformation of elastic bodies under the application of 
generalized forces. The deformation of the elastic body is computed from 
the strain energy of the body. It is an energy based approach and therein 
lies its simplicity. 

Castigliano' s theorem, developed by Alberto Castigliano in 1879, is 
a method by which one can determine the deflection of an elastic body at 
the point of application of a force. If forces F A and F B are exerted on 
an elastic body at two different points, A and B, there are four 
associated deflections [Ref 2] 

(2.1) 

6 A (2-2) 

* a a^ha^a (2-3) 

^BB~^BB^B (2.4) 

where the first subscript implies the point of interest, and the second 
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implies the force of influence. The jf's are constants, and are known as 
influence coefficients. They represent the deflection of one point 
relative to the other. For example, represents the deflection of point 
A relative to point B. These influence coefficients are properties of the 
elastic member. 

Maxwell's law of reciprocity states [Ref. 2, p. 637] 



B~^BA 



(2.5) 



which implies that the deflection at point A due to a unit force applied 
at point B is equal to the deflection at point B due to the application of 
a unit force at A. 

Castigliano developed a method where the deflection due to multiple 
forces acting on an elastic body is obtained as the summation of 
deflections due to forces applied sequentially, one at a time. The final 
result is a set of equations like those found in (2.1) through (2.4). One 
must then superimpose these equations to obtain a series of equations for 
the displacements, or 8's, at the different points where the forces act on 
the elastic body. For two forces F A and F B we obtain 

< ■ 2 • 6 > 



^B^BA*f>BB=^BA F A*^BB^B 



(2.7) 



Now one needs to determine the work done by each force at each point 
of the elastic body where it acts. The work done by F A at point A is 






=— F A 

2 *A°AA 



( 2 . 8 ) 



which, after substituting equation (2.1) becomes 
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"a 



(2.9) 



Likewise, the work done by F B at point B is 

(2 - 10) 

If force F a continues acting on the body while F B is gradually 
applied, we see that there is additional work done on point A due to F B , 
namely 

"aB^AB • 11 ) 

Similarly, if force F B continues acting on the body as F A is slowly 
applied, one can see that there is additional work done on point B due to 
F a , namely 

Hba=Fb6ba (2.12) 

After substituting equations (2.2) and (2.3) into equations (2.11) and 
(2.12), and by invoking Maxwell's law of reciprocity, equation (2.5), one 
obtains 



<2.13) 

(2.14) 

The total work done on the body, or its total strain energy, if F A is 
applied before F B is the summation of equations (2.9), (2.13), and (2.10): 

fAA?l + 2 + > ( 2 • 1 5 > 
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If force F b is applied first, and then F A , the order of addition becomes 
equations (2.10), (2.14), and (2.9) 



W=C/=-i ( f a ^i*2f AB F a F a *f M Fl) 



(2.16) 



Thus, the total work done on the body is irrespective of whether F A or F B 
is applied first. 

From Equation (2.16) one finds that the displacement at point A is 
equal to 



Castigliano ' s theorem states that for any force F i acting on an elastic 
body, the deformation or deflection at the point of application of the 
force F x is 



in the direction of F if where U is the total strain energy of the elastic 
body under the application of forces. 

B. SHAPE MEMORY ALLOYS 

To cause the deformation of a helical spring one needs to apply an 
external force. Shape Memory Alloy (SMA) active elements were chosen to 
provide the necessary external forces. An SMA active element has low mass 
and a very high force to mass ratio; this attractive feature allows the 




(2.17) 



and the displacement at point B is equal to 




(2.18) 




(2.19) 
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miniaturization of the whole structure. A brief introduction of SMA's is 
imperative to have a good understanding of how the spring will be 
deformed . 

A Shape Memory Alloy (SMA) , is a metallic alloy that is given a 
certain predetermined shape at a high temperature. Once the alloy is 
cooled, it can be deformed, and will remain deformed until heated. Once 
heated above a certain temperature, the alloy "remembers" its undeformed 
shape and returns to it . 

There are many alloys that exhibit this shape memory effect. Among 
them are Ni-Ti, Ni-Ti-Cu, Cu-Al-Ni-Mn. The alloy is first shaped into its 
desired "undeformed" shape at a high temperature, when the microstructure 
is in its austenite phase. These "undeformed" shapes can vary greatly, 
but the primary shapes we are considering are those of a thin wire of a 
given length, or a rod with a given circular curve. Once formed, the 
alloy is quenched to allow the microstructure to come into its martensite 
phase. It is now ready to be deformed. 

The alloy can now be deformed by stretching it, bending it, or 
reshaping it by any one of a number of means. It will stay deformed from 
its original shape until it is once again heated up back into the 
austenite region, where it will return to its original shape. It is 
beyond the scope of this thesis to present the microscopic analysis of 
this transformation. It is merely intended to explain what an SMA is and 
how it will be used as an actuator for spring deformation. 

The next two sections provide two spring models that describe the 
positioning of the SMA actuators relative to the spring. The SMA elements 
are in their deformed states initially, and they revert back to their 
undeformed states once heated above a certain temperature. During this 
process the spring is deformed. The undeformed shape or the memory shape 
can be given to the alloy by annealing for some time at a fixed 
temperature and then by rapid cooling back to room temperature. In the 
discussion to follow it will be assumed that the SMA has already been 
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given its memory shape and attention will be focused on the design of the 
helical spring for achieving the design goal. 



C. MODEL Is HELICAL SPRING UNDER ECCENTRIC COMPRESSION 



This section considers a spring model consisting of a helical spring 
along with its actuators such that the spring can provide two rotational 
degrees of freedom besides a single degree of freedom for linear 
translation. In this model, shown in Figure 2.1, the helical spring is 
fitted with two end caps. Three SMA wire actuators are attached to the 
end caps just outside of the spring. All of the SMA wires are to be of 
the same length, spaced 120° apart. During the process in which the SMA 
wires are placed, the spring is given a small initial bias compression. 
This keeps the SMA wires taut and eliminates any slack in the wires. 




As current is applied to heat one of the SMA's, the active element 
shrinks back to its original "unde formed" length. During this process the 
other two SMA wires remain in their deformed configuration, and the top 
plane of the spring bends over by virtue of eccentric compression. 
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Fj , Fj f and F s are forces exerted on the spring by the SMA wires, and 
F 2 , F 4 , and F 6 are dummy forces. The deflection of the points of 
application of these forces can be readily obtained using Cast igliano ' s 
theorem. Position vectors, r 1 's, from an arbitrary point A to the points 
where the forces are applied are constructed. Angle 0 is a measurement 
taken from the point where F l is applied, around in a counter-clockwise 
manner. Using these position vectors, moments and torsions due to forces 
F l through F 6 are summed up at A. 

R is the radius of the spring, and L is the length of the spring. E 
in the modulus of elasticity of the spring material, and G is the shear 
modulus of that material. I is the area moment of inertia of the cross 
section of the spring coil, and J is the polar moment of inertia. The 
number of spring coils is n. 

Using these values, the total stain energy of the spring can be 
calculated. Once this has been done, Cast igliano ' s theorem is invoked, 
and the displacement of the spring at any one of the six points of 
application of the forces can be found. Knowing the relative 
displacements of the different points on the spring coil, the angle of 
deflection can be computed. 

For the spring design problem, the angle of deflection is 
predetermined. When a helical spring is chosen, R, I, J, E, G, and n are 
known. From these quantities the force required to deflect the spring is 
calculated. If this force is one that the SMA wire can exert on the 
spring, the spring has been properly designed. If not, some of the 
parameters of the spring, geometry or material, must be changed and the 
forces recomputed. When the force required to deflect the spring matches 
the force the SMA can exert, the design problem is completed. 

D. MODEL 2: HELICAL SPRING UNDER BENDING 

In the second spring model, the assumption is that the spring bends 
under the action of SMA rods. Three SMA rods are attached to the spring 
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internally along its length and placed 120o apart. Initially the SMA rods 
are in their deformed shape. When one of the SMA rods is heated it 
regains its undeformed shape and bends the spring in the process. 




The SMA rod applies a force to each of the spring coils. Each of 
these forces are assumed to have two components, one normal to the coil 
directed towards the center of curvature of the spring denoted F, and the 
other in the tangential direction denoted f. From Figure 2.2, 0 is an 
angular measurement internal to the spring coil measured from the outer 
most point on the spring, moving in a counter-clockwise direction. R is 
the radius of the spring. The radius of curvature measured to the center 
of the spring is P, and p is the radius of curvature to the inside point 
of each of the coils as they are bent. Angle <t> is the angle of curvature 
of half of the spring. L is the length of the spring, n is the number of 
coils, E is the modulus of elasticity of the spring material, G is the 
shear modulus, I is the area moment of inertia, and J is the polar moment 
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of inertia. For an arbitrary point A on the spring, r denotes the 
position vectors from A to the points of application of the forces, and a 
is the angle subtended by each coil at the center of curvature of the 
spring . 

Consider now only the top half of the spring, since the spring is 
symmetric and the bottom half is identical to the top. When an SMA rod is 
bent through the application of current, each coil of the spring is acted 
upon by two forces, F and f. Using these forces, and the position vectors 
from point A to their point of application, a summation of all of the 
bending moments and torsions at A due to the action of these forces can be 
obtained. The total strain energy due to these moments and torsions is 
calculated, and Cast igliano ' s theorem is applied to compute the 
displacements at each point of application of the forces. 

To reiterate, for the helical spring design problem, the spring's 
displacement is a given. Knowing the equations for the strain energy, and 
the material and geometry of the spring, one can work backwards to find 
the force required to bend the spring in this manner. One needs only to 
iterate using the spring material and geometry to find a force 
commensurate with the given displacement. 
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III. DEFORMATION ANALYSIS OF THE FIRST HELICAL SPRING MODEL 



Recalling the description of the first spring model described in the 
preliminaries, one finds that F, , F J# and r 5 are forces exerted by the SMA 
wires on the spring. F 2 , F 4 , and F^ are dummy forces on the spring needed 
to find the displacements at their points of application. Point A is an 
arbitrary point on the spring, and 9 is the angular measurement from F L 
around in a counter-clockwise direction. 




Figure 3.1 Helical spring under eccentric compression. 



There are six position vectors from the six points on the spring where 
the forces are applied to point A. 



r l =R ( 1 -cos0 ) i-Rs in 0 j 



(3.1) 



r 2 = -i?(cos0-cos6O° ) i+R ( sin6 0° -sin0) j 



(3.2) 



f 3 = -R(cos6 0° *cos0) _T+i? (sin6 0° -sin0) J 



(3.3) 




r 4 = -i?(l+cos0) l-Rsin&j 



(3.4) 



f 5 = -i? (cos60° +cos0) l-R (sin60° +sin0) j 



(3.5) 



.? 6 = -\R (cos0-cos6O° ) i-i? (sin0+sin6O° ) j 



Summation of moments about point A yields 

E + g ( ^i> =0 



(3.6) 



(3.7) 



where , 

F^-Fii c i=l, 2, . . . , 6 (3.8) 

and where M* and My are the x and y components of the reaction moments at 
point A. 

After computing the cross products of ( r x x FJ , one can set the like 
vector components equal to each other and find 

M -R [- [F 1 +F 4 ) sin0+ (F 2 +FJ (sin60° -sin©) 

14 23 ( 3 . 9) 

- {F 5 +F 6 ) (sin60° +sin0) ] 



M y -R [ -F x ( 1-COS0) +F 4 ( 1 +COS0) +(F 2 +F 6 ) (COS0-COS6O 0 ) 
+ (F 3 +F S ) (cos0+cos6O°) ] 



(3.10) 



To compute the strain energy, one needs to find the bending moment and 
torsion, m and t, on the spring coil. These can be obtained from the 
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reaction moments M* and My through a coordinate transformation in the 
following manner: 




where 



r(6) = 



cos0 sin0 
-sin© cos0 



Following the coordinate conversions, the moment and tor*. 



become 



(-^+^+(-^ 2 +^+^-^) cos 60°) sin0 
+ (F 2 +F 3 -F 5 -F 6 ) sin6O°cos0] z n 



f=i?[ (f 1 +f 2 +f 3 +f 4 +f 5 +f 6 ) + (-f 2 -f 3 +f 5 +f 6 ) sin6O°sin0 
+ (-F 1 +F 4 + (-F 2 +F 3 +F 5 ~F 6 ) cos60°) cos0] z t 



The shear force at any point of the spring coil can be 
the static equilibrium of forces: 

V s (f 1 +f 2 +f 3 +f 4 +f 5 +f 6 ) £ 

The expression for strain energy due to shear is 



(3.11) 



(3.12) 



ion equations 



(3.13) 



(3.14) 



obtained from 



(3.15) 
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(3.16) 




fv 2 

2GA 



Rdd 



where U v is the energy due to shear, f is the form factor (10/9 for a solid 
circle), V is the shear force, L is the length of the spring (2nn ) , G is 
the shear modulus, A is the cross-sectional area of the spring wire, and 
x is the integration variable denoting the length of the spring coil. 
After making the substitution of J, the polar moment of inertia, for the 
circular area A by the expression 



A= 



2J 



I 



2 



and integrating with respect to 0 from 0 to 2n7l, the energy due to shear 
becomes 



U v ~- 



10rmr 2 V 2 R 

18 GxJ 



(3.17) 



One can now see that since r is much smaller that R, the r 2 term will 
dominate the numerator and make U v very small. Equation (3.17) will be 
compared later with the strain energy due to that of torsion to show this 
difference . 

The strain energy due to the bending moment is obtained through the 
equation 




(3.18) 



where U m is the strain energy due to the bending moment, m is the bending 
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moment, L is the length of the spring (2n7T) , and El is the bending 
stiffness of the spring. After integration one obtains 




(3.19) 



The strain energy due to torsion is obtained through the equation 



1 ik R< * 

0 0 



L ^ 2rm 

r t-z . r 



(3.20) 



where U t is the strain energy due to torsion, t is the torsion, L is the 
length of the spring, and JG is the torsional stiffness of the spring. 
Through integration the strain energy due to torsion becomes 



A comparison of the strain energy due to shear and the strain energy 
due to the torsional moment U t demonstrates how much smaller the strain 
energy due to shear is. Assume that their are eight coils, a spring 
radius of 5 mm, a spring coil radius of 0.5 mm, that forces F x through F 6 
are applied uniformly with a unit magnitude of 1 N, and that stainless 
steel is used as the spring's material. The ratio of strain energy due to 
torsion, U t , to the strain energy due to shear, U v , is 



u t = ( 3 F 3 2 + 5 F x F 2 + 3 F x F 3 + 2 F x F 4 + 3 F x F 5 + 5 F x F s + 3 F 2 2 
+5F 2 F 3 +3F 2 F 4 +2F 2 F s +3F 2 F # +3F 3 2 +5F 3 F 4 +3F 3 F s 
+ 2 FjF s +3F 4 +5F 4 F S +3F 4 F 6 +3 F s 2 +5F s F 6 +3 F* 2 ) 



(3.21) 
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(3.22) 



9 /*V2 (3F*+5F 1 F 2 + ~+3Fj) _ lQn 
U v 10\ zj ( F 1 +F 2 + — +F$) 2 

Since this ratio is so large, one can see that strain energy due to shear 
is very much smaller than the strain energy due to torsion. When the 
spring deforms by bending, the same can be shown for the ratio of the 
strain energy due to bending to the strain energy due to shear. This 
implies that the strain energy due to shear can be neglected. 

The total strain energy is simply the addition of equations (3.19) and 
(3.21) , 



u= U m +U t =22gl [ +F , f 2 -F, F 3 -2F x F 4 -F 3 F 5 *F x F 6 +F 2 2 +F 2 F 3 

-F 2 F 4 -2 F 2 F 5 -2F 2 F 6 +Fj + F 3 F 4 _ FjF 5 -2F 3 F s 
+ Fd + F 4 F S -F 4 F 6 +F 5 2 + F 5 F e +Fg 2 ] 

+ [3F 2 + 5F 1 F 2 +3F 1 F 3 +2F 1 F 4 +3F 1 F 5 

+ 5 F x F t + 3 F 2 + 5 F 2 F 3 + 3 F 2 F 4 + 2 F 2 F s + 3 F 2 F s 
+3F 2 +5F 3 F 4 +3F 3 F 5 +2F 3 F s +3 F 4 
+ 5 F 4 F s + 3 F 4 Fg + 3 Fl * 5 F s F 6 + 3 F 6 2 ] 



By invoking Cast igliano ' s theorem the six displacements at the points 
of application of the forces are obtained as: 



17 



(3.24) 



6l = i fir [2 Fj +F 2 -F 3 -2F 4 -F s +F s ] 

+ [6i ! i+5F 2 +3F,+2F 4 +3F 5 +5F < ] 



= [F 1+ 2F 2+ F 3 -F 4 -2F 5 -2F s ] 

(3.25) 

+ i fj^ f 5F i +6 F 2 + 5F 3+ 3F 4+ 2F s +3F s ] 



t- F i +F 2 +2F 3 + F4-F 5 -2F s ] 

(3.26) 

+ i ^ [3F 1 + 5F 2 +6F 3 +5F 4 +3 F s + 2F s ] 



= t- 2F i- F 2 + F 3 +2F 4+ F 5 -F 6 ] 

(3.27) 

+ t 2F i +3F 2 + 5F 3+ 6F 4+ 5F s+ 3F < ] 



6s= 'W f- F i- 2F 2-F 3 +F 4+ 2F 5+ F s ] 

(3.28) 

+ [ 3F i +2F 2 + 3F 3+ 5F 4+ 6F s+ 5F 4 ] 
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(3.29) 



t, * = I j£r lJ? i ~2F t -2F 3 -F A *F s *2F t ] 

+ t 5i? i + 3F 2 +2F 3 +3F 4+ 5F 5+ 6F 6 ] 



As stated before, F 2 , F 4 , and F 6 are dummy forces, which means that 
they do not really exist, they are used merely to compute the 
displacements at those points. Setting F 2 = F 4 = F 6 = 0 one obtains: 

[2F i- F 3~ F s) + -f^ [6Fi+3F 3 +3F s ] (3.30) 

b ^ B i§T [F 1 +F 3-2F S ] + [5F 1+ 5F 3 +2F 5 ] (3.31) 

*3 = J f£r [ - F 3 +2F 3 ~ F s ] [3F 1+ 6F j+ 3F 5 ] (3.32) 



= [ -2Fj +F 3 +F s ] [2Fj +5F 3 +5F S ] 

ZJG 



2 El 



(3 .33) 



= [-F 1 -F 3 +2F s ] + [3F 1 +3F 3 +6F 5 ] 

ZJG 



2 EI 



(3.34) 



a«= nU -- [Fj-2F 3 +F 5 ] + [5F,+2F 3 +5F S ] 

ZJG 



2 EI 



(3.35) 



Since all of the SMA actuators will be identical, the force they will 
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exert is expected to be the same. Consequently, if all of the three 
actuators are activated simultaneously, the deformation at the different 
points on the spring will be the same. This implies that the spring will 
undergo compression. The interest here lies in the bending of the top 
plane of the spring. Therefore, all of the three SMA actuators will not 
be activated simultaneously. Now define P as the angle by which the top 
plane of the spring bends through when one or two of the SMA wires are 
activated simultaneously. 




The two possible cases to cause the bending of the top plane of the 
spring are the application of one force, and the application of two 
forces. In the first case one force, F w is applied. This single force 
will cause displacements to occur at all six points on the spring. Of 
interest are the deflections at points 1 and 4 to determine the bending 
angle of the top plane of the spring. The deflections at points 2, 3, 5, 
and 6 will be such that the bending will not take place at in any other 
direction. By knowing the displacements at points 1 and 4, one can 
calculate p. With F L applied one finds 
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(3.36) 



From Figure 3 



Figure 3.3 Angi 



which means 



When any one 



5 4 =n7ii? 3 F 3 



f— ) 

■\ EI JG) 



(3.37) 



.3 one can see thatsimple geometry one can see that 




tanP = 



2R 



(3.33) 



P=tan 







(3.39) 



force, say F l# is applied on the spring, the angle (3 is 



given by 






(3.40) 



The other important case is when two forces are applied 
simultaneously, say F x and F 3 . In this case the important displacements 
are at points 2 and 5, where 



fi 2 =rmi? 3 {F x +F 3 ) 



\ 2EI 2 JG) 



(3.41) 



fi 5 =rmi? 3 (F x +F 3 ) 



-1 
2 El 



3 

2JG 



(3.42) 



This implies that for the application of the two forces F 3 and F 3 






(3.43) 



Since each of the SMA wire actuators will be identical, the force that 
each will exert on the spring will be equal in magnitude. This means, 
then, if two SMA wires are actuated the resultant angle of deflection will 
be 



p 



=tan‘ 1 



TmR 2 F[— +— ] 
\ El JGj 



(3.44) 



From these results one finds that if one SMA wire is actuated, or two 
are, the angle of deflection is the same. The difference is the direction 
in which the helical spring bends over. It is now clear that by 
activating one or two actuators it will be possible to bend the top coil 
of the spring in three separate directions by a positive or a negative 
angle (}. The case of actuating all three wires has not been investigated 
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since this would produce the same result as that of applying a single 
axial force through the center of the spring, or applying a normal load. 

Now that the effects of helical spring geometry and material 
properties are known in computing the deflection angle (3, use of these 
relationships will be instrumental in designing the proper helical spring. 
SMA wires deform by a constant amount of almost 5% of their original 
length in the presence or absence of external forces. However, for this 
property of the SMA to be exhibited repeatedly, the opposing stresses in 
the SMA wires should not exceed a certain value. The design of the 
helical spring will be computed from a known value of p that will give the 
force required to produce that deflection. If these forces produce 
stresses that are below the stress limit of the SMA wires, then the design 
is feasible. However, one should try to increase the stresses in the wire 
as high as possible, within limits, such that the spring to be designed is 
not too soft. A spring that is too soft will not be strong enough to act 
as a manipulator for surgical applications. 
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IV. DEFORMATION ANALYSIS OF THE SECOND HELICAL SPRING MODEL 



The second helical spring model considered is one where the spring is 
bent over directly by an SMA rod. Three rods will actually be placed 
inside the spring 120° apart to give the spring a full range of motion. 
Here, consideration is given to only one SMA rod for the development of 
the model. 




The only major assumption made is that the length of the spring is the 
same in the bent over configuration as in the undeformed configuration. 
With this in mind, certain quantities are defined. Let R be the radius of 
the spring and L be the length of the undeformed spring. In the deformed 
configuration, it is desired that the spring bends into the shape of a 
circular arc. Let P be the radius of curvature of the central axis of the 
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spring, and p the radius of curvature of the inner part of the spring. 
Let O be the angle between the top coil of the spring and the central 
coil, <(> be the angular change of spring coils, 9 be the angular 
measurement along the length of the coils of the spring, a be the angle 
subtended by two adjacent spring coils, and n be the total number of 
coils. From these definitions, one sees from Figure 4.1 that the 
following identities hold: 



* = 



L 

2P 



(4.1) 



a = 



20 

n 



(4.2) 




(4.3) 



p=P -R (4.4) 

In this model there are three different coordinate frames that are 
used. The traditional 1, ], and k frame is fixed on the top coil of the 
spring. The ±' , , and k' frame is relative to each succeeding coil. 

The £ n and E t frame is a normal and tangential transformation of the 1, 
and k frame. 

Suppose that as the SMA bends over, it exerts two forces on each 
spring coil it touches, F in the normal direction, and f in the tangential 
direction. The line of action of the normal forces passes through the 
center of curvature. For the ease of computation, it is helpful to 
translate the point of application of the normal force along its line of 
action to the center of curvature of the spring. Point A, at the top 
center of the spring, is the point where position vectors from all other 
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points on the spring will be constructed to. In considering only the 
first coil, there are two position vectors from each point on that coil, 
one corresponding to the normal force F, and one to the tangential force 
f. These position vectors have been defined this way so that they will 
apply to each coil for these forces. They are 

f 1 - ( P-RcosQ ) l '-RsinOj ' (4.5) 

z 1 '= {P-RcosQ) l ' -Rsin&j -p_T (4.6) 

The primes on the i and j vectors indicate that they are in the spring 
coil's frame of reference, where the unprimed vectors are in the fixed 
frame of the top coil. These vectors will give the position from any 
point on the first coil relative to point A. After recognizing that 

F 1 = Fjf = FiCOsfl'-FiSinQJc' (4.7) 

f 1 = '+f 1 cos4>£ ' (4.8) 

the total moment about point A caused by these forces can be calculated by 

M Ai *(g l x? l )+[Z l 'x? l ) (4.9) 

For the second coil the position vectors remain the same, but the 
forces acting on this coil are F : , f 1# F 2 , and f 2 where 

F 1 =F 1 l=F 1 cos (a+$) f '-F a sin (a-nj)) £' (4.10) 

f 1 =F 1 £=f 1 sin(a+<J>) l'+^cos (a+<f>)£' (4.11) 
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F 2 =F 2 cos$i ' -F 2 sin$>R ' (4.12) 

F 2 =f 2 sin4)2 ' + f 2 cos<t>£' (4.13) 

Now the total moment about point A is 

A f A = (f^Fj + (f x X x/j) + (i^xFg) + (f x 'x* 2 > (4.14) 

By proceeding in a like manner, one arrives at the k-th coil, where 
this coil is acted upon by the forces F 1# F 2 , . .., F k , and f 1# f 2 , ..., f k . 

These forces take the form 

F l =F l cos [ (k-1) a+4>] l'-F^ in[ (k- 1) a+<fj] R' (4 . 15) 

F 2 =F 2 cos [ (k-2) a+<()] l '-F 2 sin [ (k-2) a+4>] R' (4.16) 

F^Fj^cosil)! '-F^sinifc# ' (4.17) 

[ (k-1) a +<()] i'+^cos [ (k-1) a+<J>] R' (4.18) 

sin [ (k-2 ) a+4>] i ' + f 2 cos t (£-2) a + 4>] R ' (4.19) 

£ k =f k sin$i '+f k cos$R' (4.20) 

Now the total moment about point A becomes 
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(4.21) 



k 



k 




(f x 'xF ± ) 



or, after computing the cross products 



k 




+ [(P-J*cos0) (F i sin [ (ic-i) a+4>] 

-f ± cos [ (ic-i) a+4>] ) +f i p] j ' 

+J?sin0 (i^cos [ (ic-i) a+4>] +ijSin [ (ic-i) a+4>] ) 



The easiest way to work with this moment is to break it into 
follows 



Jc 




+f ± (JteinOcos [ (ic-i) a+4>] ) } 



M y ' = y^ iFi ( - (P-.Rcos0) sin [ (ic-i) a+4>] ) 






+f r i ( (P -RsinQ) cos [ (ic-i) a+4>] -p) ) 



k 

M z *=^2 {f 1 (-i?sin0cos [ (ic-i) a+$] ) 



+ ^ ( -i?sin0sin [ (ic-i) a+4>] ) ) 



(4.22) 



parts as 



(4.23) 



(4.24) 



(4.25) 
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A coordinate transformation is now in order to transform the moments 



from the coil fixed frame to a normal-tangent ial frame. The 

transformation matrix T(9), from equation (3.16), will be used such that 

cos0 sin0|/ x '\ (4.26) 

-sin0 cos0 \y I 




The moment equations (4.23) through (4.25) become 

k 

M n =^ iFi ( -PsinOsin [ (ic-i) a+4>] ) 

i=l 

+f d (PsinOcos [ (ic-i) a+<J)] -psin0) } 




(F 1 ( (P-Pcos0) sin t (ic-i) a+4>] ) 



+fj ( - (P-Pcos0) cos [ (ic-i) a+4>] -psin0) } 



(4.28) 



Jc 




+f i ( -PsinOsin [ (ic-i) a+<J>] ) ) 



With each of these moments in the proper frame of reference, the 
strain energy due to each of them must be calculated. With the subscript 
k on each of the energies to remind one that this is the strain energy of 
the k-th coil, the strain energy due to each of these moments is given by 
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( 4 . 30 ) 



2ic 




0 



Mn 
2 El 



RdB 



2 * 



C/ c ,*=J -z-^RdO 



2JG 



2 l m 2 

o’.-.,-/ ** 



2 El 



( 4 . 31 ) 



( 4 . 32 ) 



Since each of the moment terms are expressed as a sum of k quantities, 
the integrals involving the square of the moments becomes quite 
cumbersome. To simplify the computation, one needs to use matrix algebra. 
First define the incremental force vector as follows 
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Then define 



/A iTj\ 
AF, 



af= 



AF k 

Afi 

Af, 



e a 2 



( 4 . 33 ) 



A £ 



*■/ 






r -Psin0 [ (ic-1) o+4>] 1 

-PsinOsin [ ( k-2 ) a+<J>] 



-Psin0sin4> 



e a* 



( 4 . 34 ) 



PsinOcos [ (£-1) a+d>] -psin0^ 
PsinOcos ( (k-2) a+<p] -psin0 



Psin0cos4>-psin0 



e R* 



( 4 . 35 ) 



It can now be shown that 
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Likewise for the torsional moment, define 



f (R- PcosB) sin [ (k-1) a +4] > 
(J?-Pcos0) sin [ lk-2) a +4>] 



(£-Pcos0) sin4> 



r - (i?-Pcos0) cos t (ic-1) a +6] -pcos©^ 
- (£-Pcos0) cos [ lk-2) a+4>] -pcos0 



- (jR-Pcos0) cos<J)-pcos0 ; 



Then, 






M t =F T 



\ B » 



(4.39) 



For the moment in the z' direction, define 
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(4.40) 



r -i?sin0cos [ (k- 1) cc+<t>] ^ 
-i?sin0cos [ (k-2) a+0] 



^ -/?sin0cos<J> > 



f -Rsin6 sin [ (ic-1) a+d>] > 
-RsinQsin f (k-2 ) a +<J>] 



(4.41) 



-i?sin0sin<{> 



Then, 



M 2 *=F : 



in 



V C 2 



(4.42) 



With the moments defined by Equations (4.36), (4.39), and (4.42), one 
can obtain their squares as follows: 



M„ 2 =AF r - 



V 

A 



-AF 



(4.43) 



M 2 =AF 



T. 






A/ 



•(Bj :B 2 ) -AF 



(4.44) 
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•(CjICj) -A F 



(4.45) 



M*-=AF r - 



/ C. x 



V C 2 / 



Now define 



/a \ 



4 = 



l*a. 



•(AiiAj) = 



A A 5 AA 

AA 5 AjAj 



•A.i : A.2 
•Ajl S A 22 



one arrives at 



M„=AF t -A-AF 



where the elements of A n are of the form 



A 11 (i,j)=P 2 sin 2 0 sitl [ {k-i) a+4>] sin [ (k-j) a+4>] 



the elements of A 12 are of the form 



A 12 (i / j)=-P sin 2 0 sin [ (/c-i) a+4>] cos [ (ic-j) a+<|>] 
+Ppsin 2 0 sin [ (ic-i) a+4>] 



Ai (i.j) =A 12 (i, j) t 



and the elements of A 72 are of the form 



(4.46) 



(4.47) 



(4.48) 



(4.49) 



(4.50) 
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^22 (i, j) =P 2 sin 2 0cos [ ( k-i ) a+4>] cos [ (k-j) a+4>] 
-Ppsin 2 0(cos [ (k-i) a+4>] +cos [ ( k-j ) a+<t>] ) 

+ p 2 sin 2 0 



Similarly, 



by defining 



B= 



'B\ 



\ B 2) 







i b,b 2 




'Bu 


b 12 


•(Bi:B 2 ) = 


b 2 b. 


5 b 2 b 2 


— 


b 21 


: B 22 



one can see that 

Mc=AF t -B-AF 



where elements of B u take the form 

B ri (i, j) = (J?-Pcos0) 2 sin [ (k-i) a+4>] sin [ (k-j) a-Nj>] 



the elements of B 12 take the form 

B 12 (i t j) =- (1? - Pcos0) 2 sin [ (k-i) a+4>] cos [ (k-j) a+<J>] 
- (£-Pcos0) pcosOsin [ (k-i) a+<J>] 



B 21 (i,j)=B 12 (i,j) r 

and the elements of B 22 take the form 



( 4 . 51 ) 



( 4 . 52 ) 



( 4 . 53 ) 



( 4 . 54 ) 



( 4 . 55 ) 



( 4 . 56 ) 
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(4.57) 



B 22 (i, j) = (R- Pcos0) 2 cos [ {k-i) a+4>] cos [ {k-j) a+$] 
+ (j?-Pcos0) p cos0 (cos [ {k-i) a+4>] +cos [ {k-j) a+$] ) 

+ p 2 cos 2 0 



By defining 



c= 


( c ] 


* ( Cl : c 2 ) “ 


qq 






^11 5 ^12 




\ C 2j 


p2 C \ 


; c 2 q 




p22 ! ^22. 



(4.58) 



one can show that 

MZ'=AF T -C-AF (4.59) 

where elements of C u take the form 

C 1:1 (i, j) =R 2 sin 2 0cos [ (k-i) a+<J>] cos [ (k-j) a+4>] (4.60) 

the elements of C 12 take the form 

q 2 (i» j) =J? 2 sin 2 6sin [ (k-i) o+<J>] cos [ (k-j) a+4>] (4.61) 

C 2i d’j) =c iz O-'j) T (4.62) 



and the elements of C 22 take the form 

C 22 J) =-R 2 sin 2 0sin [ {k-i) a+4>] sin [ {k-j) a+<J>] (4.63) 



Now that the squared moments are in compact form, the integrals can 
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be evaluated for the computation of the strain energies, 
energy in the k-th coil due to the bending moment M„ is 



The strain 



2 x 




0 



m; 



2 El 



?-Rcfi = AF T - 



2 n 



dQ 



■A F 



( 4 . 64 ) 



where one must evaluate the integral of A term by term. This means that 
if one defines 




2k 

2 n 

— fFs in 2 0sin [ (ic-i) <x+4>] sin [ (k-j) a+4>] c® 

2 El j 
o 

sinacos [a (2ic-i-j+l) ] + — cos [a (i-j) ] 

a(a 2 -47i 2 ) 2 



( 4 . 65 ) 



A, 2 'u l j)uf^y 12 a,j ) dQ 



■kfr) 



2* 3 F 



a (a 2 -4n 2 ) 

«F 



sino sin [o (2k-i-j +1) ] 



( 4 . 66 ) 



sin [a (i-j) ] 



a ( 



327t 3 Pp gin/J.jginfa (jc-i+i]} 
o 2 -16it 2 ) \ 2 / 2 
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An '(i.J) aa 12 '(i. j) r 



( 4 . 67 ) 



=(-^r) { — ; T ' T - 27 sing cos t« (2k-i-j*l) ] 

\ 2FX/ a(a 2 -47i 2 ) 

+-^cos to(i-j) ] + np 2 

+ — 6 *rc 3p P --sin/-£)cos t-f (2Jr-i-j+l) J cos [-£ (i-j) ]) 
a(a 2 -l6Ti 2 ) \ 2 / 2 2 



then, the energy in the k-th coil due to the moment M„ is 



U ntk =bF T A ' *A F= A F r * 




*AF 



The energy in the k-th coil due to the torsional moment 

form 



( 4 . 68 ) 



( 4 . 69 ) 



takes the 
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(4.70) 



2it 

r 

"■•‘■/t 5 



2JG 



RdQ^AF 



T. 






dO 



•AF 



where a term by term integration, as before, yields Equation (4.75), where 
the elements of B u ', B 12 ', B 21 ' , and B 22 ' are defined as follows: 



' (i ' j) 1 2^) <7I (* 2+ t) cos to u ~ j) ] 

-f Ml - Isina cos [a (2k-l-j+l) ] ) 

1 - ’ ’ a (a 2 -4ji 2 ) J 



(4.71) 



« 2 -n 2 



Bl2 ' (i ' J,)= (2^) <7l ( i?2 + T) sinta(i ' : '' )1 

2^ +iE F(«i^lt ina sin [a (2k-l-j+l) ] 

V a a 2 -7i 2 a(a 2 -47i 2 ) ) 

J 471 (tt 2 -87i 2 )Pp _ 47tgi?p ) ain /iL\ sin [o (Jc-i + i) ] } 

^ a (o 2 -167i 2 ) a 2 -4n 2 J \2/ 2 



(4.72) 



B 21 ' (i.j) =B 12 ' (i,j) T 



(4.73) 



and 
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( 4 . 74 ) 



B 22 'U,j) =( 2?c) {7t ( i?2 + ir) cos ta (i_J) 1 +7tp2 

S -. 2 7 lg ^ + . I!^ i S _ jJ?^ . > lsinttcQS [a (2k-i-j+l) ] 
V a a 2 -7i 2 (a 2 -47i 2 ) ) 



{ 



471 (tt 2 -87I 2 )Pp _ 47ttti?p 
a(a 2 -167i 2 ) a 2 -47i 2 

47t (tt 2 -87i 2 )Pp _ 47ittigp 
a (a 2 -167t 2 ) a 2 -47t 2 



jsin|-|Jcos [a (k-i+|) ] 
jsin|-|Jcos [a (ic-j+|) ] } 



Thus the strain energy of the k-th coil due to the torsional moment M t is 



U Ct *=AF t *.B'AF=AF 



T. 



*11 ' : B 12 ' 



*21 ' : * 22 ' 



A F 



( 4 . 75 ) 



The strain energy of the k-th coil due to the moment M z , is 



U z'.k= 



2k 

f 



2EI 



Rdd=AF T ■ 






Cdd 



af 



( 4 . 76 ) 



where the elements of the submatricies C n ', C 12 ', C 21 ', and C 22 ' are defined 
in the following manner: 
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sina cos [a (2/c-i-j+l) ] 



( 4 . 77 ) 






& { 2ti 3 

2 El) a (a 2 -47i 2 ) 



+ -^-cos t® (i - J) 1 ) 



c 12 'U.j) = 



<- 



271 3 



2£J/ a ( a 2 — 47t 2 ) 



sina sin [o (2/c-i-j+l) ] 



( 4 . 78 ) 



+ ^-sin [o (i-j) ] } 



C 21 ' (i, j) =C 12 ' (i, j) r 



( 4 . 79 ) 



and 



'(i. j)=f— 



2 71 2 



V 2 ^/ a(a 2 -47i 2 ) 



sina cos [a (2ic-i-j + l) ] 



( 4 . 80 ) 



+ —■ cos [a (i-j'J ] ) 



Thus the energy of the k-th coil due to the moment M z , is 



41 



(4.81) 



U z - k =AF T -C'-AF=AF T - 



c ' \ c 

* °12 



^21 : ^22 



•AF 



The total energy of the k-th coil of the helical spring is the 
summation of the energies due to the three moments, or 

u k= u n.k+ u t.k +u z'.krbF T -(A'+B'+C') A F (4.82) 



If one defines 



DlA'+B'+C' 



(4.83) 



the incremental energy of the helical spring up through the k-th coil now 
becomes 



k k 

A (AF t 'D'AF)j (4.84) 

Cast igliano ' s theorem is now invoked to determine the incremental 
displacements of each of the coils in the normal and tangential 
directions . 
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(4.85) 




=D 11 AF 1 +2D 12 AF 2 +2D 13 AF 2 + . . . +2D lk AF k 



M 



Ac; 
A F, 



)' 



2D 21 AF 1 +D 22 AF z +2D 23 AF 2 + . 



• + 2 D 2k AF k 



(4.86) 






=2D kl AF 1 +2D k2 AF 2 +2D k3 AF 2 + . 



* D uA F k 



(4.87) 



which simplifies to 



f 8 ll 




■Du 


2D 13 


2 D 23 


• 2D lk 


'af 3 






- 


2D 2i 


&22 


2 D 23 


■ 2 D 2k 


Af 2 


(4.88) 


A, 




' CO 

£ - 
M 


CO 

£ fa 


CO 

£ - 

u> 


• Ac* . 


A F * t 





Now define 





■Djj 




2£„ • 


• • 2i? JJt 


A & 


2D 2i 


D 22 


2X>„ • 


• • 2 D 2k 




’ CO 
£ 


2 Ac* 


2-Djci • 


• • A*. 



then the final important relationship becomes 
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Af=A _1 6 



( 4 . 90 ) 



With this relationship, knowing the physical geometry and material 
properties of the helical spring, and specif ing a desired displacement, 
one can determine the force required by the SMA active element to bend the 
spring. The matrix A contains a conglomeration of information based on 
all of these known quantities. Once these given quantities have been 
specified, it is a simple matter to determine the required force matrix F. 
If the force required cannot be achieved by the SMA, iteration is required 
for the spring's geometry or material properties. A program is enclosed 
that will aid in the design of a helical spring under these conditions. 
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V. 



SUMMARY AND RECOMMENDATIONS 



This thesis investigated the deformation of helical springs under 
eccentric loading using two different models. The goal was to develop an 
algorithm for computing the deflection of a helical spring under a known 
force. The desired final shape of the spring was that of a spring bent 
over in a circular fashion. The two models used achieved this by two 
different methods, both utilizing SMA active elements in two different 
configurations . 

The development of this study was to obtain relationships between 
force and displacement so that one could predict or control the 
deformation of the spring. A spring under normal loading is governed by 
an equation such as 

F=k6 (5.1) 

where force F is proportional to displacement 5. This proportionality 
constant, k, is based on the spring geometry and material properties. In 
the development of this study, similar relationships were obtained, but 
with a more complex form of proportionality factor. This factor is 
dependent on the spring configuration, and is not a constant. The 
proportionality factor is also a function of the geometry and material 
properties . 

This study has been a first step in the design of a robotic 
manipulator to be composed of a series of helical springs that utilize SMA 
actuators to control the manipulator. In the course of future research, 
SMA actuator design is needed, and in assembling the helical spring with 
the actuators. Once the actuators are designed, the results acquire by 
this study can be employed to match a spring's geometry and material to 
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the force that can be generated by the SMA to bend the spring in a desired 
manner . 

When the design process is completed, a helical spring and actuator 
system will be realized, and can be implemented in a minimally invasive 
surgical manipulator. The conclusions of this study can also be used in 
any other situation where an eccentric load is placed on a helical spring 
causing it to bend over. 
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PROGRAMS AND PLOTS 



The first algorithm delevoped was for the first model, where the 
spring was deformed by applying an eccentric load axially on the spring. 
It is a FORTRAN program where material properties, spring geometry, and 
desired deflection angle are required as inputs, and the force required to 
deform the spring is the output. 



PROGRAM THESIS 

★ 

* This program is for the spring design used in the manipulator skelton 

* for minimally invasive surgical applications. Given the physical 

* dimensions and material properties of the spring, and the desired 

* deflection where the force will be applied, it will solve for 

* the deflection of the spring on the opposite side of where the force 

* is applied, the force required, and for the angle of deflection. 

* 

* Inputs: Bending stiffness of coil El, torsional stiffness of the coil 

* JG, number of coils n, radius of the spring R, radius of the spring 

* coil rc, length of spring, and deflection d where the force will be 

* applied. 

* Outputs: deflections d, deflection angle B, and force F. 

* 

REAL F, E, I, J, G, R, rc , L, pi, c, El, JG, dF, dOPPF, B, Z 
INTEGER n, m, k 
CHARACTER* 20 MATL 

★ 

PRINT* 

PRINT*, ' THIS PROGRAM WILL AID IN THE DESIGN OF A SPRING FOR USE' 
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PRINT*, ' IN THE SKELTON FOR MINIMALLY INVASIVE SURGICAL' 



PRINT*, ' APPLICATIONS. ' 



pi = 3.1415926 



5 PRINT* 

PRINT*, 'SELECT AND ENTER THE TYPE OF MATERIAL YOU WISH TO USE' 
PRINT*,' (i.e. STAINLESS STEEL): ' 

READ* , MATL 
PRINT* 

PRINT*, 'ENTER MODULUS OF ELASTICITY (E) , IN ENGLISH UNITS' 
PRINT*, ' (i.e. 28 . 0e6 ) : ' 

READ*, E 

PRINT* ,' ENTER SHEAR MODULUS (G) , IN ENGLISH UNITS' 

PRINT*, ' (i.e. 10 .6e6) : ' 

READ*, G 
PRINT* 

PRINT*, 'ENTER “1“ IF YOU HAVE DIMENSIONS IN ENGLISH UNITS (in), ' 
PRINT*,' OR "2" IF YOU HAVE DIMENSIONS IN SI UNITS (min): ' 

READ* , k 

PRINT*, 'ENTER RADIUS OF SPRING: ' 

READ*, R 

PRINT* ,' ENTER RADIUS OF SPRING COIL: ' 

READ*, rc 

PRINT*, 'ENTER LENGTH OF SPRING: ' 

READ* , L 

PRINT*, 'ENTER DESIRED DEFLECTION: ' 

READ*, dF 

10 IF (k.EQ.2) THEN 

R = R* (0 . 001) *39 .37 
rc = rc* (0 .001) *39 .37 
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L = L* (0.001) *39.37 
dF = dF* (0 .001) *39.37 
ENDIF 
PRINT* 

PRINT* , ' ENTER NUMBER OF SPRING COILS: ' 

READ* , n 
PRINT* 

PRINT*, 'ENTER 'T IF ONE FORCE IS APPLIED,' 

PRINT*, ' OR M 2 M IF TWO FORCES ARE APPLIED: ' 

READ*, m 

I = 0 .25*pi* (rc**4) 

J = 0 . 5*pi* ( rc**4 ) 
c = n*pi* (R**3 ) 

El = E* I 
JG = J*G 
IF (m.EQ.l) THEN 
F = dF/ (c* (1/ (El) + 3/ (JG) ) ) 
dOPPF = c*F*(-l/(EI) + 1/ ( JG) ) 

B = ( AT AN ( ( dF - dOPPF) / ( 2 *R) ) ) * ( 180 /pi ) 

ELSE 

F = dF/ (2*c* (1/ (2*EI) + 5/ (2*JG) ) ) 
dOPPF = 2*c*F*(-l/ (2*EI) + 3/ (2*JG) ) 

B = ( AT AN ( ( dF - dOPPF) / (2*R) ) ) * (180/pi) 

ENDIF 

PRINT* 

PRINT*, 'FOR MATERIAL ' , MATL 
PRINT* 

PRINT*, 'DEFLECTION AT POINT OF APPLICATION (in): ' , dF 

dF = dF/ (0 .001*39 .37) 
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PRINT*, 'DEFLECTION AT POINT OF APPLICATION (mm): ' , dF 

PRINT* 

PRINT*, 'DEFLECTION AT POINT OPPOSITE APPLICATION (in): ' , dOPPF 

dOPPF = dOPPF/ (0.001*39 .37) 

PRINT*, 'DEFLECTION AT POINT OPPOSITE APPLICATION (mm): ' , dOPFF 

PRINT* 

PRINT*, 'DEFLECTION ANGLE (degrees) : ' ,B 

PRINT*, 'FORCE NECESSARY (lbs) : ' ,F 

PRINT* 

PRINT*, 'IF YOU WISH TO BEGIN AGAIN, TYPE ANY NUMBER AND HIT' 
PRINT*,' RETURN. ELSE ENTER "99“ TO QUIT. ' 

READ* , Z 

IF (Z.NE.99) GO TO 5 
END 



The second algorithm developed was for the second model, that of a 
spring bent over directly by an SMA active element. It is a MATLAB 
program where one must input the geometry and material properties of the 
helical spring, as well as the final desired deflection angle of the 
spring. It generates a proportionality matrix whose size is based on the 
number of coils of the spring, and computes the force required to bend the 
spring. It plots force versus angle deflection to obtain a relationship 
between the magnitude of the force required and the amount of deflection. 

% thesis. m 
% 

% This program is for the spring design used in the manipulator skeleton 
% for minimally invasive surgical applications. Given the physical 
% dimensions and material properties of the spring, and the desired 
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% deflection of the spring, it will solve for the force required to bend 
% the spring using an SMA rod inside the spring. 

% 

% Inputs: Spring length (L) , radius (R) , coil radius (r) , Modulus of 

% Elasticity (E) , Shear Modulus (G) , number of coils (n) , 

% desired angle of deflection (phi). 

% Outputs: Force required to bend spring (F) . 

% 

clear 

% 

L = input (' Length of spring in millimeters '); 

R = input (' Radius of spring in millimeters '); 
r = input (' Radius of spring coil in millimeters '); 

E = input ( 'Modulus of elasticity in GPa '); 

G = input ('Shear modulus in GPa '); 

n = input ( 'Number of spring coils (must be even number) '); 
phifinal = input ('Final desired deflection angle in degrees '); 

% 

L = L/1000; 

R = R/1000; 
r = r/1000; 

E = E * 1 e 9 ; 

G = G*le9; 

phifinal = phi f inal *pi/ 180 ; 

% 

I = 0.25*pi*r^4; 

J = 0 .5*pi*r A 4; 
mom = R/ ( 2 *E* I ) ; 
tor = R/ (2*J*G) ; 

% 

for i = l:n 
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in it f ( i ) = 0 ; 
end 

f = initf ' ; 

% 

count = 0 ; 
dradius = 0; 
radius (1) = 50; 
for z = 1:200 

anglephi = L/ ( 2*radius ( z ) ) ; 

if anglephi <= phi final, count = count + 1; end 
dradius = -0 . 05*radius ( z) ; 
radius (z+1) = radius (z) + dradius; 
end 
% 

dP = 0; 

P (1) = 50; 
for q = 1: count 
% 

for i = l:n 
for j = l:n 
D(i,j) = 0; 
end 
end 

% 

ro = P (q) -R; 

phi (q) = L/(2*P(q) ) ; 

a = (2*phi (q) ) /n; 

cl = (2*pi~3)/(a* (a^2-4*pi^2) ) ; 

c2 = pi/2; 

c3 = P(q)^2; 

c4 = ( 8*pi^2 ) / (a /s 2-4*pi^2 ) ; 
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c5 = ( 64*pi"3 ) / (a* (a"2-16*pi"2) ) ; 
c6 = pi* (R"2 + (P(q) "2/2) ) ; 

c7 = ( ( R"2 *pi ) /a - ( 2*pi *a*R*P (q) ) / (a"2-pi"2 ) ... 

+ pi* (a"2-2*pi"2) *P (q) "2/ (a* (a"2-4*pi"2) ) ) ; 
c8 = ( (4*pi* (a"2-8*pi"2) *P(q) *ro) / (a* (a"2-16*pi"2 ) ) 
- ( 4 *pi *a*R*ro) / (a"2-4*pi"2) ) ; 
c9 = (32*pi"3) / (a* (a"2-16*pi"2) ) ; 
clO = sin (a) ; 
ell = sin (a/2 ) ; 
cl2 = mom*R"2 ; 
cl3 = cl*c3*cl0; 
cl4 = c2*c3; 
cl5 = P (q) *ro*c9*cll ; 
cl6 = c7*cl0; 
cl7 = cl*clO; 
c!8 = P (q) *ro*c5 ; 



for k = l:n/2 
for j = l:k 
for i = 1 : k 



cl9 = cos (a* (2*k-i- j+1) ) ; 
c20 = sin (a* (2 *k-i - j +1 ) ) ; 
c21 = cos (a* ( i- j ) ) ; 
c22 = sin (a* ( i- j ) ) ; 
c23 = s in (a* (k-i+0 . 5 ) ) ; 
c24 = sin (a* (k-j+0 . 5) ) ; 



All ( i , j ) = mom* (cl3 *cl9 + cl4*c21); 

A12(i,j) = mom*(cl3*c20 - cl4*c22 - cl5*c23); 
A21(i,j) = mom*(c!3*c20 + c!4*c22 - c!5*c24); 
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A2 2 ( i , j ) 



% 



% 



% 



= mom* ( -cl3 *cl9 + cl4*c21 + pi*ro^2... 

+ cl8*cll*cos ( (a/2)*(2*k-i-j+l) ) *cos ((a/2)*(i-j))) 



Bll ( i , j ) 
B12(i, j) 
B21 ( i , j ) 
B22 ( i , j ) 



= tor* ( c6*c2 1 - c 1 6 * c 1 9 ) ; 

= tor*(c6*c22 - cl6*c20 + c8*cll*c23); 

= tor* ( -c6*c22 - cl6*c20 + c8*cll*c24) ; 

= tor* ( c6*c2 1 + cl6*cl9 - c8*cll* (cos (a* (k-i+0 . 5 ) ) 
+ cos (a* (k- j + 0 . 5) ) ) + pi*ro /s 2); 



Cll (i, j ) 
C12 ( i , j ) 
C21 ( i , j ) 
C22 ( i , j ) 



cl2* (c2*c21-cl7*cl9) ; 
cl2* (c2*c22-cl7*c20) ; 
cl2* (-c2*c22-cl7*c20) ; 
cl2* (c2*c21+cl7*cl9) ; 



end 

end 



dll = All+Bl 1+Cl 1 ; 
dl2 = A12+B12+C12 ; 
d21 = A21+B21+C21 ; 
d22 = A22+B22+C22; 
for i = l:k 
for j = l:k 

D ( i , j ) = D ( i , j ) + dll ( i , j ) ; 

D(i+k,j) = D(i+k,j) + d21(i,j); 

D ( i , j + k) = D ( i , j +k) + dl2(i,j); 

D(i+k,j+k) = D(i+k,j+k) + d2 2 ( i , j ) ; 
end 
end 
end 

D = 2 . * D ; 
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for i 



1 : n 



D(i,i) = 0 . 5 *D ( i , i ) ; 
end 

% 

for i = l:n/2 

theta(i) = (n/2 - (i-l))*a; 
deltan(i) = (1 - cos ( theta ( i ))) *dP ; 
deltat(i) = ( sin ( theta ( i ) ) - theta ( i )) *dP; 
end 

% 

delta = [deltan' ;deltat ' ] ; 

F = inv(D) *delta; 
for i = 1 : n/2 

fx(i) = F ( i ) *cos ( theta ( i ) ) - F (n/2+i ) *sin ( theta ( i ) ) 
fy(i) = F ( i ) *sin ( theta ( i ) ) + F ( n/2+i ) *cos ( theta ( i ) ) 
end 

fxy = [ fx ' ; f y ' ] ; 
f = f + fxy; 
force ( : , q) = f ; 

PHI (q) = phi (q) * (180/pi) ; 
dP = -0 . 05*P (q) ; 

P(q+1) = P (q) + dP; 
end 
% 

clg 

for j = 1 : n/2 

plot ( PHI , force ( j , : ) ) 

title ( 'X-Force vs Radius Angle Phi, for Coil #') 
xlabel ( ' Phi (degrees ) ' ) 
ylabel ( 'Force (N) ' ) 
meta thesplot 
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end 



for j = n/2+l:n 

plot ( PHI , force ( j , : ) ) 

title ( 'Y-Force vs Radius Angle Phi, for Coil #') 
xlabel ( ' Phi (degrees ) ' ) 
ylabel ( 'Force (N) ' ) 
meta thesplot 
end 



The final element is a set of plots using the previous program. They 
are force versus deflection angle graphs generated by the program using 
the inputs of a 20 mm long spring, with a radius of 5 mm, coil radius of 
0.5 mm, and eight coils. The material chosen was stainless steel with a 
modulus of elasticity of 190 GPa, and shear modulus of 73 GPa. The final 
desired deflection angle was 30o . The eight plots are for the top half of 
the spring. The forces in the bottom half are going to be identical as 
the spring is symmetrical. Of the eight force/deformation plots, four are 
for normal forces, and four are for tangential forces. 
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